Mixed neutron-star-plus-wormhole systems: Linear stability analysis 
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We consider configurations consisting of a neutron star with a wormhole at the core. The wormhole 
is held open by a ghost scalar field with a quartic coupling. The neutron matter is described by a 
perfect fluid with a polytropic equation of state. We obtain static regular solutions for these systems. 
A stability analysis, however, shows that they are unstable with respect to linear perturbations. 

PACS numbers: 04.40.Dg, 04.40.-b, 97.10. Cv 



O ' I. INTRODUCTION 

(N 

' The geometrical model of electric charge suggested by Wheeler in the middle of the 1950s, describes a tunnel 
^ . connecting two space-time regions that is filled by an electric field This idea has stimulated many studies of 
solutions with a nontrivial space-time topology - wormholes. One of the most significant contributions in this area 
is the model of a traversable Lorentzian wormhole, suggested by Morris and Thorne Q • The traversability assumes 
04 ■ that matter and radiation can travel freely through the wormhole. The key condition for the existence of such a type 
of wormholes is the necessity to violate the weak/null energy conditions. In Einstein gravity this means, that the 
• matter creating the wormhole must possess very exotic properties. 

Thus the question arises whether such exotic forms of matter can exist in the Universe. A strong argument in favor 
, of such a possibility is the observed accelerated expansion of the present Universe. The cause of this acceleration 
bJQ' is attributed to a repulsive component of the Universe, the so-called dark energy. Its amount exceeds 70% of the 
total energy density of the Universe. With its large negative pressure, dark energy causes the Universe to expand 
■ increasingly fast, according to a power law or exponentially. Moreover, astronomical observations (see, e.g., Refs. [HHI 
| and more recent estimates Q ) indicate the possibility, that an even more exotic form of energy exists in the Universe, 
f — ■ called phantom dark energy. The presence of such an energy assumes the violation of the weak/null energy conditions 
and results in even faster acceleration. 

If dark energy, in one form or another, does indeed exist in the Universe, this would give the basis for concluding that 
localized compact objects consisting of dark energy might also exist. Such objects could be so-called dark energy stars 
CN . [||-[l5j - objects consisting of dark energy and possessing a trivial topology. On the other hand, if phantom dark energy 
is present in the Universe then it is possible to imagine a situation where topologically nontrivial configurations - 
wormholes, can exist p^ - l2lj . The possibility is not excluded that the sizes of such configurations could be quite large 
and comparable to those of various astrophysical objects - such as ordinary stars or neutron stars p2j . 

One type of matter providing violation of the weak/null energy conditions is the so-called ghost or phantom scalar 
field. With the opposite sign in front of its kinetic energy term, such a field allows for solutions with nontrivial 
topology, including traversable wormholes. The properties of such solutions as, for instance, their regularity and 
stability then depend on the particular field employed. Perhaps the simplest possibility to obtain wormhole-like 
solutions is to consider massless scalar fields [HI, [24[. However, as recently shown, such wormhole configurations are 
unstable with respect to linear [2^, [26| and nonlinear perturbations [l?} . 

When stepping beyond the bounds of Einstein gravity and considering the issue of the stability of wormhole solutions 
within the general framework of scalar-tensor theories with massless non-minimally coupled scalar fields it is seen that 
such solutions are also unstable under linear perturbations [28l - l3lj ]. 

In order to obtain stable solutions, the study of configurations supported by a ghost scalar field with a self-interaction 
potential might seem promising. For wormholes such potentials were considered in Refs. [32l. l33j. In particular, here 
a ghost scalar field with a quartic coupling was employed, and it was concluded that regular, stable solutions with 
topologically nontrivial (wormhole-like) geometry exist. A similar conclusion was reached in Ref. [34| for wormholes 
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with a sine- Gordon ghost scalar field. 

Here our objective is the study of configurations with non-trivial topolog y c onsisting both of ordinary and exotic 
matter. We have considered such mixed gravitating systems before in Refs. JH, where we suggested the possible 
existence of configurations consisting of a traversable wormhole (supported by a massless ghost scalar field) filled by 
a perfect polytropic fluid. We have shown that static, regular solutions can indeed be constructed, which describe 
such mixed star-plus- wormhole systems. These possess new physical properties which distinguish them from ordinary 
stars. 



In [35| we made some preliminary estimates concerning the stability of such mixed configurations with respect to 
linear perturbations. However, our analysis was incomplete since it was performed only in the external region of the 
star. On the other hand, bearing in mind that a wormhole supported by a massless scalar field without ordinary 
matter is unstable j25l [26j , we may naively expect that adding ordinary matter to such a configuration will not lead 
to the stabilization of the system. One reason could be that the main contribution to the e nerg y density near the 
throat is coming from the scalar field, as our studies of mixed configurations performed in [3^, indicated. 

As mentioned above, in [2^, [2(| the question of stability of wormhole-like solutions with a massless ghost scalar field 
was clarified. In particular, it was shown that when allowing for perturbations of the throat radius the solutions are 
unstable with respect to spherically symmetric perturbations. Recent investigations in (37j showed that wormhole-like 
solutions remain unstable for special choices of the scalar field potential. 

Bearing all this in mind, in the present paper we reanalyze the stability of wormhole solutions for a ghost scalar 
field with a quartic coupling, first addressed in Ref. (33j . Unlike those pioneering calculations, we now allow for per- 
turbations of the throat radius. Subsequently, we perform such a linear stability analysis for the mixed configurations, 
consisting of a wormhole and both ordinary matter and a ghost scalar field with a quartic coupling. 

The paper is organized as follows: In Sec. Ill Al the general set of equations is derived for configurations consisting of 
a neutron star with a wormhole at the core, where the neutron matter is modeled by a perfect fluid with a polytropic 
equation of state. In Sec. IIIB1 we present numerically obtained static solutions for such topologicallly non-trivial 
configurations supported by a ghost scalar field with a quartic coupling. In Sec. IIIH a linear stability analysis is 
performed for these solutions. Finally, in Sec. IIVI our results are summarized. 

II. STATIC CONFIGURATIONS 
A. General equations 

We consider a model of a gravitating ghost scalar field in the presence of a perfect fluid. The Lagrangian for this 
system is chosen as 



L = --^—R-\d ll tpd^-V(<p)+L rn . (1) 
107rG I 



Here <p is the ghost scalar field with the potential V(<p), and L m is the Lagrangian of the perfect isotropic fluid (where 
isotropic means that the radial and the tangential pressure of the fluid agree) , which has the form L, n = p [3a . l39j . 
Using this Lagrangian, the corresponding energy-momentum tensor can be presented as 



T t k = (e + P ) Ul u k - b k lV - di(pd k ip - 8$ 



(2) 



where e and p are the energy density and the pressure of the fluid, u l is the four-velocity. The metric can be taken in 
the general form 

ds 2 = e v {dx°f - e x dr 2 - e»dQ?, (3) 

where v, A, and fj, are functions of the radial coordinate r and the time coordinate x° = ct, and dfl 2 is the metric on 
the unit two-sphere. 

In considering equilibrium wormhole-like configurations, it is convenient to use the polar Gaussian coordinates 

ds 2 = e"{dx°f - dr 2 - e^dn 2 , (4) 
where now v and \l are functions of r only. Introducing the new function R defined by e M = R 2 , the (§), (1), and (|) 
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components of the Einstein equations with metric (j4]) take the form 
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where the prime denotes differentiation with respect to r. 

The equation for the scalar field ip resulting from the Lagrangian (fl} is 
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Using the metric (d}, this equation gives 
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Not all of the Einstein field equations are independent because of the conservation of energy and momentum, T^ k = 0. 
Taking the i — 1 component of this equation gives 



R' 



0. 



(10) 



Taking into account the expressions 



Tn = T? 



3 = -P 



and also Eq. ([9]), we obtain from Eq. ([TO 
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(11) 



To model the matter filling the wormhole, it is necessary to choose an appropriate equation of state. In doing 
so, we proceed from the assumption that our configuration is essentially a relativistic object, where the wormhole is 
filled with relativistic matter having a pressure comparable with its energy density. For this kind of matter we choose 
neutron matter. In much of the literature neutron matter is described by more or less conventional equations of state, 
reflecting general properties of neutron matter at high densities and pressures. Various forms of such equations of 
state can be found, for instance, in Refs. j40l - |43l |. 

Since in the present paper we consider only general properties of mixed neutron-star-plus-wormhole systems, we 
restrict ourselves to a simplified variant of the equation of state, where a more or less realistic neutron matter equation 
of state is approximated in the form of a polytropic equation of state. Namely, we employ the following parametric 
relation between the pressure and the energy density of the fluid, 



P 



7-1 



kc n 



rh) 



where rib is the baryon number density, is some characteristic value of m, m b is the baryon mass, and k and 7 
are parameters whose values depend on the properties of the neutron matter. 
It is convenient to rewrite the above equation of state (EOS) in the form 



rs l+l/r 

P = Kp b 



£ = p b C 



np, 



(12) 



with the constant K = fee 2 '{n^ 11 ' ) m b ) 1 1 , the polytropic index n = 1/(7 — 1), and p b = n b m b denotes the rest-mass 
density of the neutron fluid. 



4 



Setting m& = 1.66 x 10 _24 g and — 0.1 fm~ 3 , we consider below configurations with k = 0.1 and 7 = 2 [44| . 

corresponding to a gas of baryons interacting via a vector meson field, as described by Zel'dovich [45|, [46| (see also 
Ref . [47j where relativistic configurations with such an equation of state were considered) . In Ref . [36[ we have already 
considered a similar mixed system consisting of a massless ghost scalar field and a neutron fluid with an EOS in the 
form of (|12[) . Here we extend those results to the case where the scalar field potential is present. 

Introducing the new variable 6 [46| 



Pb = PbcO n 



(13) 



where pb c is the density of the neutron fluid at the wormholc throat (or, in other words, at the core of the configuration) 
we may rewrite the pressure and the energy density, Eq. (|12[) , in the form 



PbcC 



nK Pbc 



Making use of this expression, we obtain for the internal region with 9^0 from Eq. (jlll) 

dO dv 
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where a = Kp^J 1 /c 2 = p c /(pb c c 2 ) is a constant, related to the pressure p c of the fluid at the wormhole throat. This 
equation may be integrated to give in the internal region with 6^0 the metric function e v in terms of 9, 



1 + a{n + 1) 



1 



1 



(16) 



and e v " is the value of e v at the throat where 9 = 1. The integration constant v c is fixed by requiring that the 
space-time is asymptotically flat, i.e., e u = 1 at infinity. 

Thus we have three unknown functions - R,9 and ip - for which there are four equations ([5])- ([7]) and ((9]) (only three 
of which are independent), and also the relation (|16[) . For the numerical calculations it is convenient to rewrite these 
equations in terms of dimensionless variables. Since in section lilBl wc will consider a particular case, where the scalar 
field tp is equal to zero at the wormhole throat, but its derivative is nonzero, we can introduce dimensionless variables 
as follows: The potential can be expanded in the neighborhood of the throat as 



1 3 
<Pir+ -<p 3 r , 
6 



where <pi is the derivative at the throat, the square of which corresponds to the "kinetic" energy of scalar field. Then, 
it is convenient to use new dimensionless variables expressed in units of tp\ . Namely, introducing 
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with L having dimensions of length, one can rewrite Eqs. ©-(O and 
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(18) 

(19) 
(20) 
(21) 



where V = V/ip\ is the dimensionless potential of the field, and B = (pb c c 2 ) I V? is the dimensionless ratio of the fluid 
energy density to that of the scalar field at the throat. 

Thus, the static configurations under consideration arc described by any three equations from the system (|18l) - (|2ip 
together with Eq. (fT5|) or p^|) . Note, that in the case of B = we are dealing with a system consisting of a pure 
scalar field configuration with no ordinary matter. 
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B. Quartic potential 

In this section we discuss the numerical solutions of the set of equations (p~5|) . (fT8|) - ([2T|) . We seek regular solutions 
of these equations describing configurations with a finite mass. In the case of a massless scalar field, it was shown in 
Refs. [35|, y6| that static solutions for mixed systems consisting of a scalar field and a polytropic fluid can indeed be 
obtained. 

In the present paper we consider the case where the potential has the well-known form of the </? 4 theory. It was 
shown in Refs. [321 133| that a ghost field with such a potential admits regular topologically non-trivial solutions. Our 
aim here is to study the influence which the presence of a polytropic fluid has on such solutions. 

For our purpose, we choose the potential term in the form [32l . |33j 




where m v and / are constants. Using the dimensionless variables (JTTJ) , one can rewrite this potential as follows 

V = -1/(0) (1 - A 2 2 ) 2 , (22) 

where A 2 = / 2 c 4 /(87rG) is a dimensionless constant, and V"(0) = (m^/f) 2 /(2ip1) is the value of the potential at the 
local minimum, where 4> m in = 0. Maxima of this potential are located at the points <\> max = ±1/A. Asymptotically, 
as £ — > oo, a regular solution must approach one of these maxima. Note here that, in constrast to the case of usual 
(non-ghost) scalar fields, the existence of regular solutions is only possible, when the sign of the potential is reversed 
with respect to the usual case. This is the reason that we choose a potential of the form (f22j) . that is unbounded from 
below. 

Eqs. (fT5]) - ([2"l"j) are to be solved for given a, n, and B, subject to the boundary conditions at the core of the 
configuration £ = 0, 

9(0) = 1, E(0) = E C) S'(0) = 0, i/(0) = v c , 0(0) = 0, ^(0) = 1. (23) 

The quantity S c is the eigen-parameter of the system. It is determined from the condition of obtaining asymptotically 
vacuum solutions, where <fi — > ±1/A. In turn, the value of the parameter B may be obtained by expressing ipi in 
terms of L from (jTT)) . and this gives B = 8nGpb c (L/c) 2 . Thus, the value of B is determined by the core density of the 
fluid and the choice of the characteristic size L of the configuration under consideration. The case B — > corresponds 
to the exclusion of the fluid from the system, leaving only wormholes supported by a scalar field [32l [HJ . 

Let us emphasize again that, for a given value of A, and thus a given theor y, t here is only a single value of the 
wormhole size E c for which a regular solution exists (with the potential (|2"2")) ') j32ll33l|. This situation differs from other 
variants of bosonic configurations considered in the literature. For example, in the case of systems with self- interacting 
complex scalar fields considered in [48|, |49|, the mass of the configurations, for a fixed coupling constant, is a function 
of the central value of the scalar field, and correspondingly of its central energy density. In this case one can obtain 
the dependence of the mass on the central density which, as in the case of configurations consisting only of ordinary 
(for instance, neutron) matter, may be employed in considering the stability of the configurations within the energy 
approach. However, for the wormholes of |32l . |33j such an analysis cannot be performed. 

For the mixed configurations considered in this paper, where, besides the ghost scalar field with a </> 4 potential, 
also ordinary matter is present, the total mass of the configurations depends on the amount of fluid in the system, 
and correspondingly on the fluid core density. Then the possibility of considering the stability of the configurations 
within the energy approach appears again possible and will be discussed below. 

Substituting the potential (|2"2"|) into Eqs. ([15]), ([T5[) - ([2"Tj) and using the boundary conditions (f2"3"[) , we seek a numerical 
solution of this set of equations. In doing so, the configurations under consideration can be subdivided into two regions: 
(i) the internal one, where both the scalar field and the fluid are present; (ii) the external one, where only the scalar 
field is present. Correspondingly, the solutions in the external region are obtained by using Eqs. ||18 )> -(|21 )> . in which 
8 is set to zero. The internal solutions must be matched with the external ones at the boundary of the fluid, £ = £&, 
by equating the corresponding values of the functions <fi, S, v and their derivatives. The boundary of the fluid £f, is 
defined by = 0. Knowledge of the asymptotic solutions in turn allows to determine the value of the integration 
constant v c at the throat, proceeding from the requirement of asymptotic flatness of the external solutions. 

Let us now address the total mass of the configurations. For the spherically symmetric metric (@}, the mass m(r) 
inside the radius r can be defined as follows 

r 2 Att r 

m{r) = —R c + ^J r T°R 2 dR, (24) 
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TABLE I: Characteristics of a set of configurations for EOS <|12[) . The radius of the throat R t h, the proper radius of the fluid 
flprop, and the gravitational radius r g (last column) are given in kilometers. The total mass M, the mass at the throat Mth, 
the mass of the fluid Ma, and the external part of the mass of the scalar field M s f ex t are given in solar mass units. 



Pbcg cm 3 


Rth, km R plop , km 


B 


M/Mq 


M th /M e 


M fl /M 


M sfext /M Q 


r g , km 


Without a wormhole 


l.OxlO 13 


- 


35.3940 


— 


0.2777 


- 


0.2777 


— 


0.8192 


3.0xl0 13 




34.7518 




0.7573 




0.7573 




2.2338 


5.0xl0 13 




34.1519 


z 


1.1532 




1.1532 




3.4014 


7.0xl0 13 


— 


33.5906 


— 


1.4818 




1.4818 


— 


4.3708 


l.OxlO 14 


— 


32.8113 


— 


1.8760 




1.8760 


— 


5.5335 


2.0xl0 14 


— 


30.6535 


— 


2.6482 




2.6482 


— 


7.8112 


3.0xl0 14 


— 


28.9877 


— 


2.9832 




2.9832 


— 


8.7993 


4.0xl0 14 


— 


27.6575 


— 


3.1213 




3.1213 


— 


9.2068 


6.0xl0 14 


— 


25.6722 


— 


3.1564 




3.1564 


— 


9.3101 


8.0xl0 14 


— 


24.2625 


— 


3.0779 




3.0779 


— 


9.0787 


l.OxlO 15 


— 


23.2159 


— 


2.9715 




2.9715 


— 


8.7647 


1.5xl0 15 


— 


21.5291 


— 


2.7156 




2.7156 


— 


8.0099 


2.0xl0 15 


— 


20.5768 


— 


2.5160 




2.5160 


— 


7.4212 


3.0 xlO 15 


— 


19.7367 


— 


2.2541 




2.2541 


— 


6.6487 








A = 


1,L= 10 km 








l.OxlO 13 


2.0075 


0.3165 


0.0002 


0.3823 


0.6806 


0.0000 


-0.2520 


1.1277 


5.0xl0 13 


2.0078 


0.7822 


0.0009 


0.3827 


0.6807 


0.0002 


-0.1186 


1.1288 


l.OxlO 14 


2.0084 


1.2784 


0.0019 


0.3837 


0.6809 


0.0008 


-0.0366 


1.1318 


3.0xl0 14 


2.0134 


5.5204 


0.0056 


0.4078 


0.6826 


0.0231 


-0.0000 


1.2029 


4.0xl0 14 


2.0172 


13.1029 


0.0075 


0.5194 


0.6839 


0.1335 


-0.0000 


1.5321 


5.0xl0 14 


2.0219 


21.8576 


0.0093 


0.9602 


0.6855 


0.5731 


-0.0000 


2.8322 


6.0xl0 14 


2.0276 


24.5778 


0.0112 


1.4644 


0.6874 


1.0760 


-0.0000 


4.3195 


7.0xl0 14 


2.0343 


25.0595 


0.0130 


1.8199 


0.6897 


1.4300 


-0.0000 


5.3679 


l.OxlO 15 


2.0603 


24.0895 


0.0186 


2.2767 


0.6985 


1.8818 


-0.0000 


6.7153 


1.5xl0 15 


2.1272 


22.0069 


0.0279 


2.3277 


0.7212 


1.9217 


-0.0000 


6.8658 


2.0xl0 15 


2.2318 


20.2965 


0.0373 


2.1844 


0.7566 


1.7645 


-0.0000 


6.4432 


3.0xl0 15 


2.6474 


16.4199 


0.0559 


1.8366 


0.8976 


1.3854 


-0.0000 


5.4172 


3.5xl0 15 


3.0947 


13.4552 


0.0652 


1.7418 


1.0492 


1.2971 


-0.0000 


5.1377 



where R c is the radius of the wormhole throat denned by R c = min{ R(r)}. Without loss of generality, we can take 
this throat to occur at r = 0. Note, that for the total mass the upper limit of the integral is infinity, since formally 
the energy density of the scalar field becomes equal to zero only asymptotically, as R — > oo. In practice, however, the 
scalar field decays exponentially fast. Consequently, for the values of the parameters employed all mass is concentrated 
within a size of order L. Note also, that in evaluating the above integral it is necessary to perform the calculations 
separately in the internal and external regions. 

In the dimensionless variables ([T7|) . the expression (|24l) takes the form 
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where the coefficient in front of the curly brackets has the dimension of mass. 

Using the data from Table HI the numerical results are illustrated in Figs. Q] and [2] In performing the calculations, 
we start from an initial configuration without ordinary matter (a pure wormhole) . Then, by gradually adding neutron 
matter to this system (as expressed in the growth of pb c or, equivalently, of B) we can monitor the behavior of the 
masses and sizes of the configurations. The dependence of the total mass (in solar mass units) on the core density of 
the fluid ptc (in grams per cubic centimeter) is shown in Fig. [TJ As the core density decreases, i.e., as B — > 0, the 
total mass of the configurations tends to the mass of a pure wormhole without any fluid, My/H ~ 0.38M Q . 

On the other hand, with increasing pb c , the total mass M of the configurations rises monotonically to a maximum 
and then decreases again. This is typical for this type of configurations. In Tooper's paper [47| . for instance, ordinary 
neutron stars were investigated in detail for various values of the polytropic index n. In particular, it was shown that 
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Core density of the fluid p hc , g cm" 

FIG. 1: Total mass of the configurations (in solar mass units), M/Mq, versus the core density pb c (in units of g cm -3 ) for 
EOS (|12|l with k = 0.1 and 7 = 2 (n = 1). The characteristic size L is taken as 10 km. The value of A is taken as 1. Stable 
configurations should reside to the left of the first mass peak, if they would exist. 




Total mass MIM^ 



FIG. 2: The mass-radius relation. The points are labeled with the appropriate values of the core density pbc (in units of 
). The tentative stability region is situated approximately in the range 7 x 10 14 < pb< 



10 14 g cm" 3 



< 14 x 10 14 g cm" 3 . 
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the first peak in the mass corresponds to the point dividing stable and unstable neutron-star configurations. This 
first mass peak is reached at a critical value of the core density, p^ . 

However, looking at Fig. [2j where the mass-radius relation for the configurations is presented, one can see that as the 
core density of the fluid increases, first a simultaneous growth of the size and the mass of the configurations takes place, 
which is typical for unstable compact astrophysical systems. In the interval 7 x 10 14 g cm~ 3 < pb c < 14 x 10 14 g cm~ 3 
the mass then continues to increase while the radius of the configurations decreases. Such a behavior is typical 
for stable compact astrophysical systems. Finally, for pi, c > 14 x 10 14 g cm -3 the mass and the radius decrease 
simultaneously, as it is characteristic for unstable configurations. 

Thus, the presence of the wormholc has a remarkable effect on the configurations for small central densities, since 
for ordinary neutron stars (without a wormhole) the mass increases as the radius decreases in this low density range. 
Consequently, ordinary neutron stars (without a wormhole) are stable in this range, and remain stable all the way 
to the maximum of the mass. Only beyond the associated critical value p^T a star would become unstable (see, 
e.g., [47||). The naive analogy with ordinary compact stars would then suggest, that for star-plus-wormhole systems 
instability occurs not only for pt, c > p^\ but also for small values of the core density pb c - Clearly, at this point a 
reliable stability analysis is called for. 



III. LINEAR STABILITY ANALYSIS 

The simplest configurations with a nontrivial wormhole-like topology are obtained by using a massless scalar field [23l . 
While the first stability studies (see, e.g., Ref. [llj]) did not find unstable modes, recent work showed that 
these configurations are linearly [25|, [26| and nonlinearly (2tJ unstable. In particular, the recent analyses allowed 
for perturbations of the throat, where a Schrodinger-likc equation with a regularized effective potential had to be 
solved 

In this section we perform a linear stability analysis of the above static solutions with the quartic potential (1221) . 
Also in this case it was shown, that the solutions are stable against a special type of linear perturbation, where 
oscillations of the throat radius are excluded [Hj]. Here we reanalyze the stability of these solutions, by allowing for 
perturbations that are not required to vanish at the throat. 

Subsequently, we consider the stability of star-plus-wormhole systems, where besides the ghost scalar field also a 
polytropic fluid is present. Our investigation is facilitated by the fact, that both background solutions and perturba- 
tions decay exponentially fast, when a <fi 4 potential considered [33j . 



A. General equations 



We now consider spherically symmetric perturbations of the above equilibrium configurations. In obtaining the set 
of equations for the perturbations, we will neglect all quantities which are of second and higher order. 

For the energy momentum tensor of the fluid we need the components of the four-velocity in the metric ([3]) [50| 



with the three-velocity 



—Vail vo/2 1 —va/2 Xn — vn/2 



dr 

V= d^ <<1 - 



The index on the metric functions indicates the static, zeroth order solutions of the Einstein equations. The 
components of the energy-momentum tensor @ then take the form 

T° = s-\e-»tf- l -e-\> 2 + V{<p), (26) 

Tl = -p + i e - V + \e~ V' 2 + V(<p), (27) 

Td = (e + p)uqu x - doipd 1 ^ = (e„ + p )v + e~ V if/, (28) 

Ti =T$ = -p+ V V - ±e- V 2 + V{<p). (29) 

In the above equations, the prime and dot denote differentiation with respect to r and a; , respectively. 
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Now wc consider perturbations of the static solutions. Let us denote by y any one of the functions v, A, fi, e,p or (p. 
Then we assume that y is of the form 



y = yo + y P , 



(30) 



where the index again refers to the static, zeroth order solutions, and the index p indicates the perturbation. 
Substituting these expressions into the (g), ([), and (?,) components of the Einstein equations, written in the metric 



^0 



G% 



A / 



1 



1 



A 4 + 7 m" - oM'A' + - A e~ v M + 2 A A + e 



\e- x (p> 2 + 2 M V) + e"" f/i - \iiv + ^ + 



8ttG 



— e 

4 



x (A'// - 2 M " - M ' 2 - i/// + AV - 2v" - v' 2 ) 
A/i + /i 2 - !>/i + 2/i - Az> + 2A + A 2 ) = -^T 2 2 , 



(31) 
(32) 

(33) 



we find to linear order in y p 



3 -A 



3 , , 1 



Mp + ^MoMp ~ o (Mo A p + A oM P ) ~ a p Mo + jVo ~ oMo^o 



1 >M 



e^Vp 



8ttG 



2 6 



-A 



+ ^) n' + v' q ij! p - \ p + fj,' u' 



8ttG 



Pp + e" A Vo (v£ - + V v 



1 , 1 



'(Mp + *>p) 



2 p 
16?rG 



Mp' + "p + Mo Up + o^P " o A P + M ~ 9 A P + ^ ) ~ A P 



-_Pp - e A Vo I V 



2</>o a p 



Mo +^o+o (^o 2 + «tf + M) 



(34) 



(35) 



(36) 



where V p = [d v V]o ip p - Next, from the (1-0) component of the Einstein equations 

1 



—A 



2/i' - A/u' +fi(fjf- v') 



8ttG , 



wc have to linear order 



1 

2 C 



-A 



8ttG 



2/ip' - ApMo + Mp (Mo - "a) 
Now we introduce a "Lagrangian displacement" C, with respect to x° [5 

V dx° ' 



[(eo+Po)v + e a Vp^o] 



(37) 



Then equation (1371) can be integrated to give 



\ = lh + ~j\ Mp - \v-v v 'q - -^r eX ° [i £ o +Po)C + e A V Vp] 



2 

Mo 



(38) 



In turn, the z = 1 component of the law of conservation of energy and momentum, X* fc = 0, gives 



^ v' + yl 



?i - 2 ( T 2 + T : 



0. 



(39) 
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Substituting here the components (|26l) - ([29|) . we find to linear order 



[(eo+ P o)v + e-^ p ]-^ + ^ 



-An 



(f' fj, p + 2fj,' (w' p - ^oAp^j 



The perturbed scalar field equation is found from Eq. ([8]). and is given to linear order by 

v ; - e x °- v °$ p + l(v'o-\'o + 2m' ) + \ (»p -K + H) ^ 





'd 2 V 






'dV~ 






dip 2 


ip = ip 


dip 


<P=Vo } 



= 0. 



(40) 



(41) 



Thus, we have a set of five general linear equations (|34l) . (|35|) . (|38l) . (j40|) . and (|4"Tj) for the perturbations # p , i/ p , A p , /i p , 
and iy9 p . 



B. Quartic potential 



Now we employ the set of equations obtained in the previous subsection to study the stability of the static solutions 
considered in Sec. Ill Bl In this case the potential V is given by Eq. (|22j) . and the equation of state of the fluid is 
defined by Eq. (|14|) . From these expressions one finds that the perturbed components of the pressure p p and the 
energy density e p are 



Pp = K(n + i)pl+ 1/n ei;e p , 

In turn, the static components are 



e p = np bc c 



-+a(n + l) 



y "P- 



Pa = Kpl+ 1/n e^+\ so = Pbcc 2 (1 + anOo) 0%. 



(42) 



(43) 



To proceed with the stability analysis we now assume that the harmonic perturbations have the following time 
dependence, 



(44) 



where the functions y p (£) depend only on the spatial coordinate £. For convenience, we hereafter drop the bar. 

Let us now consider the gauge freedom of the problem. First, we have the freedom of choosing the radial coordinate 
r. Here our choice has been to set Ao = 0. Second, we can make a gauge choice for the metric perturbations v p , X p , /z p . 
In particular, we may impose a relation among the metric perturbations. We are guided by the perturbed scalar field 
equation (|41j) , which contains the term (v' p — \' p + 2/i p ) ip' . This equation is considerably simplified by the gauge 
choice 



v p — A p + 2\i p = 



u p A p 2^i p . 



In the following we reformulate the general set of equations ([Ml) . (|55"j) . (|3"5|) . (|40l) and (|41l) with this choice of 
gauge, which allows us to eliminate the perturbation v p from the set of equations. Employing again the dimensionlcss 
variables dTTJ), Eq. flU) yields 



1 



4 + 5 H + 2Mo) 4>' p + u 2 e- v ^ p = -4A 2 V(0) [(1 - 3A 2 </> 2 ) ^ + (l - A 2 2 ) </> A 



Pi > 



(45) 



where we have introduced a new dimensionlcss frequency ui — ujL, and subsequently again dropped the bar for 
notational simplicity. 
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The equations (|34|) and (|36f give, respectively, 



nB 



+ <r{n + 1) 



flo «p - 0o ( ^ - 2^oA P ) + 4AV(0) (1 - A 2 2 ) O P \ , 



[Mo +^0 + ^0) 



+ w 2 e- l/ °(^ p + A p ) 



= -2 



-Btr(n + 1)0^ - ( - ^ Ap ) + 4A 2 V>(0) (l - A 2 </> 2 ) fafa 



(46) 



(47) 



while Eq. (gDJ leads to 



{stf" [1 + o-(n + l)flo] V + 0o0p} - Btr(n + l)^(^o^) + ^ 



+ 00 ( 0p - 7^0^ ) + ^0 



0p - 2 Ap + ^o A p) 



n 

-+<7(ra + l)< 

#0 



y p"0 



(48) 



-Be? [1 + a(n + 1)6 ] (a; - 2 M ;) + -^ 2 (a; - 2 M p 

0oMp + 2// [tip - \<t>'oK 



^000 (V P - T^oAp^ + 00 



= 0, 



where 



^ = 4A 2 ^(0) [(1 - 3A^ 2 ) O P + (1 - A 2 2 ) </> o 0p] • 

Note, that we can eliminate ip in (|48p by replacing the curly brackets by the corresponding expression from the 
constaint equation 



2 Mp + (/4 - !/ ) Mp - Mo A p - 2 {£?#o [1 + tr(n + 1)0 O ] ^ + O P } = 0, 



(49) 



which results from Eq. (|38|) . Here ip = C,/L is the dimensionless Lagrangian displacement. The second constraint 
equation follows from Eq. (|35p and reads 



Mo(A p - M p ) + ^o/4 - Ap Q/i 2 + Mo^oJ + 2 {oJ 2 e-»° + e"" ) Mp 



-B<r(n + + O ( ^ - -</> A p ) + 4A 2 F(0) (l - A 2 ^ 2 ) O P 



(50) 



Thus, for the four functions cf> p , X p , /i p , 9 P , we have the set of four equations (|45|) - p8|) . to investigate the stability of 
the configurations. For this set of equations, we choose the following boundary conditions at £ = 

Ap(0) = A pO , Mp(0)=Mp0, p (0) = 6 p0 , &,(0)=0, cf>' p (0) = (f> pl , (51) 
where X p , fi p , 9 p are even functions, while (p p is an odd function. The value of <p p \ is obtained from Eq. (|5T 



p i = Ba(n + l)6 p0 + -X p0 - I w e ,yc + yu p0 - 

Thus the system contains three parameters, X p o,fi p o, and 9 p $. Their values are chosen such that the following 
conditions satisfied: (i) At the boundary of the fluid, £ = the value of 9 p should remain finite to ensure that p p 
[cf. Eq. (|4"2")l ] meets the condition p p = at the boundary (see, e.g., Eq. (60) in Ref. [5(j). (ii) Asymptotically, as 
£ — > ±00, the perturbations A p , p, p , <p p should tend to zero. In this connection it is useful to determine the asymptotic 
behavior of the solutions. This can be given in analytic form: 
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FIG. 3: The case without the fluid: The lowest eigenvalue ujq is shown as a function of the throat radius E c . As E c grows, 
uj 2 -0. 



(A): Static solutions: 



90 



1/A-Ct exp 



i 



l-2C 2 /e 



8A 2 V(0)£ 
C 2 /f; 



'ft 



(B): Perturbations: 



/itp -> C 4 exp (—• \A 



A p — > — CW— w 2 £exp (— V— w 2 £). 



Here /3 2 



-8A 2 V r (0), and the are integration constants. Hence, to obtain decaying solutions for the perturbation 



(j) p of the scalar field, /3 2 should be negative. Taking into account that in our case V"(0) is always positive, this means 
that the condition uj 2 < 8A 2 y(0) should be satisfied. 

Let us now move to the results of the numerical calculations. The set of equations (|45j) - (|48|) together with the 
boundary conditions (|51j) defines an eigenvalue problem for to 2 . The question of stability is thus reduced to a study 
of the possible values of to 2 . If any of the values of uj 2 are found to be negative, then the perturbations will grow and 
the configurations in question will be unstable against radial oscillations. 

The results of the calculation of the eigenvalue u!q are shown in Figs. |3] and |U Fig.[3]shows the dependence of Wq on 
the dimensionless throat radius E c of the wormhole solutions without ordinary matter (i.e., with parameter B = 0). 
This case corresponds to the problem considered in (33|. We thus find, that the eigenvalue ujq is negative for any 
size of the throat. All these wormhole solutions are thus linearly unstable. This result is in contrast to the restricted 
previous analysis, which did not allow for perturbations of the throat radius. 

Fig. [4] shows our results for the case when the system contains the fluid, i.e., for the star-plus- wormhole systems. 
Here the eigenvalue Wq is shown as a function of the core density of the fluid p& c . As the background solutions we 
employ the static solutions obtained in section HTbI The initial value n P (0) in (|5Tj) is chosen to be /x p o = 1> an d the 
values A p o and 9 p o are chosen in such a way that the solutions exhibit the asymptotic behavior shown in (B). It is seen 
from Fig. |4]that the square of the cigcnfrcquency remains always negative, independent of p& c . Thus, unfortunately, 
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0.00 I 1 1 1 1 — i — i — i i i 1 1 1 1 — i — i — i — r-| r 




Core density of the fluid p , g cm" 



FIG. 4: The case with the fluid: The lowest eigenvalue ujo is shown as a function of the core density pb c (for L = 10 km and 
A = 1). When p bc = B = (i.e., without the fluid) co'i « -0.11610638. 

the star-plus-wormhole systems obtained with the ghost scalar field with a quartic potential are always unstable 
against linear perturbations. 

IV. CONCLUSION 

Here we continued our study of star-plus-wormhole systems, begun in [H, |3(|. The star-plus-wormhole systems 
considered there were supported by a massless scalar field. In the present paper we considered the case when the 
non-trivial topology is provided by a scalar field with a quartic potential. This choice of potential was based on the 
hope that it would allow for stable solutions. 

Choosing a neutron fluid with a poly tropic equation of state (|12[) . that is filling the wormhole, we showed that 
regular solutions for such mixed systems exist. The resulting neutron stars have a finite size, as in the case of the 
mixed configurations with a massless scalar field [35l [36j . 

For our convenience in performing the stability analysis of such star-plus-wormhole systems, we chose the values of 
the parameters in such a way that the main mass of the configurations (more than 99%) is concentrated within the 
radius corresponding to the edge of the fluid. In this case we showed that: 

(1) There exist regular static solutions found numerically by solving the coupled Einstein- matter equations subject 
to a set of appropriate boundary conditions. 

(2) Analyzing the dependence of the total mass on the core density of the fluid (see Fig. [I) and the mass-radius 
relation (see Fig. |2J), we found a range of core densities where stable configurations seemed to be possible. 

(3) The linear stability analysis performed in section UlIBl however, indicated that the square of the lowest eigen- 
frcquency of the perturbations is always negative. 

This means that all star-plus-wormhole systems considered here are unstable against linear perturbations. Moreover, 
the wormholes themselves (i.e., the solutions without ordinary matter) are unstable as well. We conclude, that the 
star-plus-wormhole systems inherit their instability from the wormholes. 

In order to find stable star-plus-wormhole systems one should thus start from stable wormholes. A possibility here 
would be to go beyond Einstein gravity and include higher curvature corrections (Hll [Hj • Such solutions would not 
need any exotic matter for their existence. 

In the static solutions considered here the ghost scalar field tends asymptotically in each of the universes to a 
different vacuum value, given by the two degenerate minima of the potential. The scalar field has therefore the shape 
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of a kink [53(, when considered as a function of the radial coordinate. A similar behavior is found in the case of the 
star-plus- wormhole configurations considered in Refs. [HI, [36|, where asymptotically the massless ghost scalar field 
assumes values that are equal in magnitude but have opposite signs. 

Another interesting possibility consists in considering wormhole solutions, which are supported by two interacting 
ghost scalar fields. An example of such solutions was given in Rcf. [54 |. Here the scalar fields assume the same value 
asymptotically in both universes. In our future work, we plan to perform a stability analysis of those solutions. 
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